Enhancing Mixing and Diffusion with Plastic Flow 
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We use numerical simulations to examine two-dimensional particle mixtures that strongly phase 
separate in equilibrium. When the system is externally driven in the presence of quenched disorder, 
plastic flow occurs in the form of meandering and strongly mixing channels. In some cases this 
can produce a fast and complete mixing of previously segregated particle species, as well as an 
enhancement of transverse diffusion even in the absence of thermal fluctuations. We map the 
mixing phase diagram as a function of external driving and quenched disorder parameters. 

PACS numbers: 05.40.-a,05.60.-k,82.70.Dd 



There have been a growing number of experiments on 
collections of small particles such as colloids moving over 
periodic or complex energy landscapes generated by var- 
ious optical methods [l], M S IE M, 0, 0] or struc- 
tured surfaces Such static and dynamical substrates 
can produce a variety of new particle segregation mech- 
anisms 0, 0, H, 0] as well as novel types of logic devices 
Driven particles on periodic substrates can also ex- 
hibit enhanced diffusive properties such as the recently 
proposed giant enhancement of the diffusion which oc- 
curs at the threshold between pinned and sliding states 

0,0, [11113, 113, [3| • This enhancement has been demon- 
strated experimentally for colloids moving over a periodic 
optical substrate and could be important for appli- 
cations which require mixing and dispersing of different 
species of particles . A limiting factor for using diffu- 
sion enhancement to mix particles is that the diffusion is 
enhanced only in the direction of the external drive. For 
instance, in a two-dimensional system with a corrugated 
potential that is tilted in the direction of the corruga- 
tion barriers, there is no enhancement of the diffusion in 
the direction transverse to the corrugation barriers at the 
pinned to sliding threshold. It would be very valuable to 
identify a substrate that allows for strong enhancement 
of the diffusion in the direction transverse to the tilt of 
the substrate, or one that would facilitate the mixing of 
particle species that are intrinsically phase separated in 
equilibrium. Such a substrate could be used to perform 
fast mixing of species and would have applications in mi- 
crofiuidics, chemical synthesis, and creation of emulsions 
and dispersions. 

In this work we show that a phase separated binary as- 
sembly of interacting particles in the presence of a two- 
dimensional random substrate tilted by a driving field 
undergoes rapid mixing and has an enhancement of the 
diffusion transverse to the tilt direction. The motion of 
the particles occurs via plastic flow in the form of mean- 
dering channels which have significant excursions in the 
direction perpendicular to the drive, leading to mixing of 
the two particle species. The mixing and diffusion occur 



even in the absence of thermal fluctuations and arise due 
to the complex multi-particle interactions. We map the 
mixing phase diagram as a function of external drive and 
substrate properties and identify regimes of rapid mixing. 
We find that as the difference between the two particle 
species increases, the mixing becomes increasingly asym- 
metric with one species penetrating more rapidly into the 
other. Our work shows that plastic flow can be used as a 
mechanism for mixing applications, and also provides a 
new system for the study of collective dynamical effects. 

We simulate a two-dimensional system with periodic 
boundary conditions in the x and y directions contain- 
ing two species of Yukawa particles labeled A and B 
with charges qa and qb, respectively. The particle- 
particle interaction potential between particles i and j 
of charges qt and qj at positions r.^ and is V{rij) = 



EoQiqj exp(- 



])lrij, where Eq = Z*^/47reeo, e is the 



dielectric constant, Z* is the unit of charge, k is the 



screening length, and 



,|. We fix K = 4/ao 



where aq is the unit of length in the simulation. The 
system size is L = 48ao. The motion of particle i is de- 
termined by integration of the overdamped equation of 
motion 



' dt 



(1) 



where rj is the damping term which is set equal to unity. 
Here F^'^ = is the particle-particle in- 

teraction force, where N is the total number of particles 
in the system. The particle density is p = N jl?. The 
substrate force Ff = — ^^(''jfe) comes from TVp 

parabolic trapping sites placed randomly throughout the 
sample. Here V^irih) = -{Fp/2rp){rik - rp)'^Q{rp - r^k), 



where Fp is the pinning strength. 



0.2ao is the pin 



radius, r.ik = jr^ — r[?'''| is the distance between particle 
i and a pin at position r[,^^ , and O is the Heaviside step 
function. The pin density is pp = Np/L/^. The external 
driving force F^; — is applied uniformly to all the 
particles. The units of force and time are Fq — Eo/ao 
and T = fl/Eo, respectively. We neglect thermal flue- 
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FIG. 1: Red circles and red lines: particle positions and 
trajectories for species A; blue circles and blue lines: particle 
positions and trajectories for species B; open black circles: 
pinning site locations in a system with particle density p = 
0.7, pin density pp = 0.34, and pinning force Fp = 1.0 at 
different driving forces, (a) Fd = 0.0, (b) Fd = 0.1, (c) Fd = 
0.4, and (d) Fd = 1.1. 



tuations so that T = 0. If the two particle species are 
initialized in a phase separated state, in the absence of 
an external drive and disorder the particles will not mix 
unless the temperature is raised above melting. 

In Fig. 1(a) we show the initial phase separated parti- 
cle configuration for a 50:50 mixture of the two particle 
species with qA/qs = 3/2 and qa ~ 3. The particles are 
placed in a triangular lattice of density p — 0.7 which 
is immediately distorted by the pinning sites of density 
Pp — 0.34 and strength Fp — 1.0. Species A occupies a 
larger fraction of the sample due to its larger charge qA 
and correspondingly larger lattice constant compared to 
species B. An external driving force Fd is applied in the 
x-direction and held at a fixed value. 

Figure 1(b) illustrates the particle trajectories at Fd ~ 
0.1 over a period of 10^ simulation steps. The trajecto- 
ries form meandering riverlike structures with significant 
displacements in the direction transverse to the drive, 
producing intersecting channels that permit species A to 
mix with species B. When the trajectories and parti- 
cle positions are followed for a longer period of time, 
the amount of mixing in the system increases. The 
riverlike channel structures are typical of plastic flow 
of particles in random disorder, where a portion of the 
particles are temporarily trapped at pinning sites while 
other particles move past, so that the particles do not 
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FIG. 2: Measure of local homogeneity H vs time for the sys- 
tem in Fig. 1 at Fd = 0.05 (black), 0.1 (red), 0.25 (green), 
0.4 (blue), 0.5 (brown), 0.6 (magenta), 0.7 (violet), and 1.1 
(top curve). H — 1 for phase segregation and H = 0.5 for 
complete mixing. 



keep their same neighbors over time. This type of plas- 
tic flow has been observed in numerous one-component 
systems including vortices in type-II superconductors 
S [3Jl3, E, Ea mi , electron flow in metal dot ar- 
rays J22] , and general fluid flow through random disorder 
2a . 24| . These works have shown that by changing the 
strength and size of the disorder, the amount of trans- 
verse wandering or tortuosity of the riverlike channels 
can be adjusted, and that these channels appear even for 
T = [ifl 



17, 1 



19| . In our system we measure the dif- 
fusion in the y-direction, dy = \{ri{t) • y — ri(0) • y}P, 
and flnd a long time transverse diffusive motion with 
dy{t) oc t" and a = 1.0, indicative of normal diffusion. 
Single component systems exhibiting plastic flow also 
show a similar transverse diffusive behavior Jl7|. The 
diffusion in our system is not induced by thermal motion 
but rather occurs due to the complex many-body parti- 
cle interactions that give rise to the meandering riverlike 
channels. In Fig. 1(c) we plot the particle trajectories 
in the same system at Fd = 0.4. At this drive, a larger 
fraction of the particles are mobile and the riverlike chan- 
nels become broader. As the drive is further increased, 
all the particles are depinned, the meandering river like 
structures are lost, and the mixing of the particles de- 
creases. Such a state is shown in Fig. 1(d) at Fd = 1.1. 
For higher values of F^; > 1.1, ffow similar to that shown 
in Fig. 1(d) appears. 

In order to quantify the mixing, for each particle we 
identify the closest neighboring particles by performing 
a Voronoi tesselation on the positions of all particles in 
the system. We then determine the probability H that a 
particle is of the same species as its neighbors. If the sys- 



3 



—i \ lJ_J \ lJ_J \ I ' ' ^ ^ ' ' \ lJ_J \ L_kj I J-l— lJ 



0.2 0.4 0.6 0.8 1 1.2 1.4 

Fp 

FIG. 3: Mixing phase diagram of pinning density pp vs pin- 
ning strength Fd in the form of a height map of the local 
homogeneity H obtained from a series of simulations with 
Fp = 1.0 and particle density p — 0.7. Strong mixing regions 
are blue and weak mixing regions are red. 

tem is thoroughly mixed, the local homogeneity H = 0.5, 
while if it is completely phase separated, H is slightly less 
than one due to the boundary between the two species. 
In Fig. 2 we plot H{t) for the system in Fig. 1 at dif- 
ferent values of Fd ranging from Fd = 0.05 to Fd = 1.1. 
For the lower drives Fd < 0.1, there are few channels 
and a portion of the particles remain pinned through- 
out the duration of the simulation so that mixing satu- 
rates near H — 0.6 to 0.7. For the intermediate drives 
0.1 < Fd < 0.5 any given particle is only intermittently 
pinned, so at long times all the particles take part in the 
motion and the system fully mixes, as indicated by the 
saturation of H to H = 0.5. For drives 0.5 < < 0.9 
the system can still completely mix but the time to reach 
full mixing increases with Fd- At Fd > 0.9 where the par- 
ticles are completely depinned, the mixing becomes very 
slow as shown by the H{t) behavior for Fd = 1.1. Within 
the strongly mixing regime, H(t) oc Aexp{—t) at early 
times before complete mixing occurs. 

In Fig. 3 we plot the mixing phase diagram of pinning 
density pp versus driving force Fd as determined by the 
local homogeneity H obtained from a series of simula- 
tions with Fp = 1.0 and p = 0.7. The value of H is 
measured after 3 x 10^ simulation time steps. Blue indi- 
cates strong mixing and red indicates weak mixing. For 
> 1.0 and all values of pp, all of the particles are mov- 
ing in a fashion similar to that illustrated in Fig. 1(d). 
Since the plastic flow is lost, mixing is very inefficient 
in this regime. For Fd < 0.6 at high pinning densities 
Pp > 0.7, most of the particles are pinned, preventing a 
significant amount of mixing from occurring. A region of 
strong mixing appears at 0.6 < Fi < 0.9 for all values of 
Pp. Here, the particles intermittently pin and depin, pro- 
ducing the large amount of plastic motion necessary to 
generate mixing. There is another strong region of mix- 




FIG. 4: (a) Black line: the average particle velocity V vs 
Fd for a system with fp — 1.0, Pp — 0.34, and p = 0.7. Red 
line: the corresponding dV/dFa curve, (b) Black circles: lo- 
cal homogeneity H; red squares: net transverse displacement 
dy for the same system as in (a). The high mixing regime 
{H < 0.6) is correlated with enhanced transverse displace- 
ments and the peak in dV/dFd- dy has been shifted down 
for presentation purposes, (c) Particle positions (circles) and 
trajectories for species A (black) and species B (blue) in a sys- 
tem with qA/qB = 3 and Fd ~ 0.2. The mixing is asymmetric 
with species A moving into the region occupied by species B 
before species B moves into the area occupied by species A. 

ing for lower pinning densities 0.2 < pp < 0.4 and low 
Fd < OA. In this regime there are more particles than 
pinning sites so that interstitial particles, which are not 
trapped by pinning sites but which experience a caging 
force from neighboring pinned particles, are present. At 
low drives the interstitial particles easily escape from the 
caging potential and move through the system; however, 
the pinned particles remain trapped so that the intersti- 
tial particles form meandering paths through the pinned 
particles. This result shows that even a moderately small 
amount of disorder combined with a small drive can gen- 
erate mixing. As the pinning density is further decreased 
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to pp < 0.15, the amount of mixing also decreases. 

In Fig. 4(a) we demonstrate how the mixing phases are 
connected to the transport properties of the system by 
plotting the net particle velocity V — {N~^ Sili ' ^) 
and dV/dFd versus driving force Fd for a system with 
Pp — 0.34 and Fp = 1.0. Here is the velocity of particle 
i. In Brownian systems, it was previously shown that an 
enhanced diffusion peak is correlated with a peak in the 
derivative of the velocity force curve [13, EH El, E3, E3] ■ 
Figure 4(a) shows that there is a peak in dV/dFd span- 
ning 0.5 < Fd < 0.9 which also corresponds to the region 
of high mixing in Fig. 3. There is also a smaller peak in 
dV/dFd at small drives Fd < 0.2 produced by the easy 
flow of interstitial particles. For Fd > 1.0, V increases 
linearly with Fd since the entire system is sliding freely. 
In Fig. 4(b) we plot the local homogeneity H for the 
same system taken from the phase diagram in Fig. 3. 
The maximum mixing {H < 0.6) falls in the same region 
of Fd where the peak in dV/dFd occurs. Figure 4(b) also 
shows that the net traverse particle displacement dy has 
peaks in the strong mixing regimes. 

We have also examined the effect of significantly in- 
creasing qA/qn so that the system is even more strongly 
phase separated. In general, we find the same mixing fea- 
tures described previously; however, the time required 
for complete mixing to occur increases with increasing 
Qa/qb- The mixing also becomes asymmetric: the more 
highly charged species A invades the region occupied by 
species B before the less highly charged species B spreads 
evenly throughout the sample. In Fig. 4(c) we illustrate 
the particle trajectories during the first 3 x 10^ simula- 
tion time steps for a system with qA/qs = 3 at F^; = 0.2. 
The mixing asymmetry can be seen from the fact that 
the black trails corresponding to the motion of species A 
overlap the blue trails representing the motion of species 
B, but the region originally occupied by species A con- 
tains no blue trails. 

One issue is whether the results reported here apply 
more generally for other types of particle interactions. 
We considered only Yukawa interactions; however, the 
meandering channel structures which lead to the mix- 
ing are a universal feature of one-component systems un- 
dergoing plastic flow though random quenched disorder. 
Studies performed on systems with long-range logarith- 
mic interactions 17| as well as short range interactions 



[24l which show this plasticity lead us to believe that 
plastic flow generated by random disorder can produce 
enhanced mixing for a wide range of particle interactions. 
For our specific system of Yukawa particles, experiments 
on single component systems have already identified a 
channel- like plastic flow regime 

In summary, we have shown that two-dimensional plas- 
tic flow induced by quenched disorder in the absence of 
thermal fluctuations can lead to efhcient mixing and en- 
hanced diffusion in phase separating systems. This mix- 
ing occurs due to the meandering of particles through 



river like flow structures. We map the general mixing 
phase diagram and find that mixing is optimized in 
regimes where the particles depin in an intermittent fash- 
ion. For higher external drives the mixing is strongly re- 
duced. These results should be general to a variety of 
systems where meandering flow channels appear. 



This work was carried out under the auspices of the 
NNSA of the U.S. DoE at LANL under Contract No. 
DE-AC52-06NA25396. 
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